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Kilometer-scale deep under-ice or -water Cherenkov neutrino detectors may detect muon and 
electron neutrinos from astrophysical sources at energies of a TeV and above. Tau neutrinos are 
also expected from these sources due to neutrino flavor oscillations in vacuum, and tau neutrinos 
are free of atmospheric background at a much lower energy than muon and electron neutrinos. 
Identification of tau neutrinos is expected to be possible above the PeV energy range through the 
"double bang" and "lollipop" signatures. We discuss another signature of tau in the PeV-EeV 
range, arising from the decay of tau leptons inside the detector to much brighter muons. 

PACS numbers: 



I. INTRODUCTION 

Kilometer-scale Cherenkov neutrino detectors now 
planned or under construction, such as IceCube at the 
South Pole and KM3NeT in the Mediterranean, are 
poised to detect high energy neutrinos from astrophysi- 
cal sources such as gamma-ray bursts (GRBs) and active 
galactic nuclei (AGN). Ultrahigh energy (UHE) cosmic 
rays with energies exceeding ~ 10 20 eV are thought to 
originate from these extraordinary sources of the highest 
energy 7-rays observed. High energy neutrinos should 
also be produced in these sources as the result of pho- 
tomeson (jyy) and/or proton- proton (pp) interactions of 
shock accelerated protons with ambient radiation fields 
and/or plasma material. 

After propagating virtually unimpeded through the 
universe, high energy neutrinos detected by experiments 
on Earth may reveal the physical and astrophysical condi- 
tions of their sources at energies and distances unmatched 
by other methods. In addition to measuring the energy 
and direction of astrophysical neutrinos, future neutrino 
telescopes will be able to distinguish between the three 
known flavors of neutrinos (and anti- neutrinos), namely 
electron, muon and tau neutrinos (y e , and v r ), by 
looking at the signature(s) of their interactions in the 
detection media. 

Tau neutrinos are not produced in appreciable num- 
bers in astrophysical sources, but will appear in numbers 
comparable to v e and as a consequence of flavor os- 
cillation between their sources and Earth. Astrophysical 
neutrino sources generically are believed to produce neu- 
trinos through 7T^~ (and K^) decay, leading to a neutrino 
flavor ratio at production of v e : : v T = 1 : 2 : 0. How- 
ever, neutrino flavor eigenstates v a (a = e, /i, t) and mass 
eigenstates Vj (j = 1, 2, 3) are mixed by a unitary matrix 
U defined as v a = V U*,Vj. The oscillation probability 
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is the phase of the slower (solar) neutrino oscillation 
and D is the distance traveled by the neutrinos before 
being detected. The oscillation probability reduces to 



>v 



E 7 Wp 3 



u, 



for ip ^> 1, which holds for 



essentially all astrophysical sources. We use the stan- 
dard expression for U a .j from Ref. [3| with solar mixing 
angle 0q = 9\2 = 32.5° and the atmospheric mixing an- 
le # a tm = 023 = 45°, following the results from the SNO 
and K2K j5| experiments, respectively. The unknown 
mixing angle 6*13 and the CP violating phase may be as- 
sumed to be zero given the current upper bounds from 
reactor experiments. The resulting astrophysical neu- 
trino flux ratio on Earth is thus generally expected to be 
v e : : v T = 1 : 1 : 1, although this may be modified 
somewhat by effects such as neutron decay, two-photon 
annihilation to muon pairs, or muon synchrotron cooling 
in the source environment 0, Q, IE @ • 

Tau neutrinos are particularly interesting because lo- 
cal backgrounds to astrophysical v T signals are low. The 
possible backgrounds are from UHE cosmic rays inter- 
acting in Earth's atmosphere and producing short-lived 
charmed mesons which deca y a s D s — > tv t , known as 
the "prompt" neutrino flux |lCt lllj . and from conven- 
tional atmospheric v e or produced in these cosmic ray 
air showers oscillating to v T as they traverse the Earth 
before being detected. 

The prompt atmospheric v T flux is not precisely known 
due to uncertainties in the extrapolation of the parton 
distribution functions to low x and in the composition 
of the high energy cosmic rays. However, the different 
models lead to expected rates of 10~ 3 to 10 -2 prompt tau 
events above 1 PeV per year in a km 3 neutrino telescope 

mm q m. 

As for oscillations to v T , the maximum propagation dis- 
tance L ks 10 4 km, the Earth's diameter, corresponds to 
a — ► v r oscillation probability V « 10 _3 ( J B I ,/TeV) -2 
(neglecting matter effects, which would reduce the oscil- 
lations), so oscillations will produce a negligible flux of 
v r for E v at the PeV scale, even in comparison to the low 
prompt fluxes. The only remaining source of PeV-scale 
tau neutrinos is extraterrestrial. 
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Tau neutrinos are also interesting because of the phe- 
nomenon of v T regeneration |16|. The high energy neu- 
trino interaction cross section rises approximately lin- 
early with energy, and the mean free path through the 
Earth becomes shorter than the Earth's diameter for 
E v > 100 TeV. The Earth is thus opaque to v e and v^ L at 
PeV energies and above, and only horizontal and down- 
going neutrinos can be observed. For v T , however, the 
produced in a charged current (CC) neutrino interac- 
tion will usually decay back to v T before losing significant 
amounts of energy, effectively regenerating the v T beam 
and leaving an upgoing v T flux up to the PeV scale. 

The best-known signature for detecting v r in a water 
or ice Cherenkov detector is called the "double bang" 
[l7j . In these events, a CC neutrino- nucleon interaction 
v T N — > tX produces a hadronic shower (denoted X), 
with the subsequent decay of the r lepton producing a 
second shower, connected to the first by the r lepton 
track. The second shower may also be hadronic, or it 
may be electromagnetic in the case of r — > ev T i7 e . The 
r produced in the CC interaction has energy (E T ) » 
§.75E V [l^, and the two showers are separated by the 
tau decay length l T = -fct T ~ 50 (E T /PeV) m (neglecting 
energy losses along the track) . Due to the short r lifetime 
and wide spacing of the detection elements in kilometer- 
scale neutrino telescopes, this signature is only expected 
to be detectable for v T with energy E v > PeV. Above 
~ 20 PeV, the typical decay length exceeds 1 km, so both 
showers usually will not be contained in a kilometer-scale 
detector; the resulting signature of a tau track and one 
shower is known as a "lollipop" |l9j . 

In this paper, we point out another distinctive signa- 
ture of extremely high energy (EHE, PeV-EeV) tau lep- 
tons, produced by the muonic decay of a r inside the 
instrumented detection volume. Although the muon has 
lower energy than the parent tau lepton, it will emit more 
light than the tau. The lepton track will thus appear 
to suddenly increase in brightness by an amount which 
should be detectable in a neutrino telescope. 

The energy range over which this signature is observ- 
able is constrained at the lower end by the requirement 
that a reasonably long tau lepton track be observed prior 
to the tau decay. At the higher end, the rising rate of 
tau photonuclear energy loss causes the brightness of the 
tau to approach that of the daughter muon above EeV 
energies. It should be noted that these limits apply to 
the energy of the r lepton in the detector; events from 
higher energy v T could be observed if the initial neutrino 
interaction vertex is some distance from the detector so 
that the r lepton loses energy in stochastic interactions 
before decaying within the detector. 



II. SIGNATURE OF MUONIC DECAY 

Identification of tau events through muonic decay, 
t — > /i^ T f M , requires the decay to occur within the de- 
tector so that the increase in brightness will be observed. 



The branching ratio for this decay channel is measured to 
be Tfj — 17.36% 3], so only a fraction of tau leptons will 
manifest themselves via this signature. However, at en- 
ergies > 20 PeV, the tau track length L T becomes longer 
than the geometric scale of the detector and double bangs 
are no longer visible. Only the lollipop and muonic de- 
cay signatures can be used to identify taus in this regime, 
even though only a fraction ~ (1 km)/L T of the taus can 
be tagged by these methods; lollipops are generally con- 
sidered to be identifiable only when the final, not the 
initial, bang is observed (see Section llV|) . 

Tau leptons are produced in neutrino V — A interac- 
tions, which at the energies of interest produce polarized 
taus. The spectrum of muon energies from the decay of 
polarized r is dn/dx = |(1 — x 3 ), where x = E^/E T 
|13tl2fil21j . The expected muon energy is thus 

= 0.4 E T . 

Depolarization via r interactions prior to decay should be 
small, and would reduce the muon energy only slightly, 
to 0.35-E T in the limit of complete depolarization |2fJ. 

Although the muon has less than half the energy of the 
tau, it appears brighter because the muon loses energy 
more rapidly than the tau, as discussed below. 



A. Energy Loss Rates 

The average energy loss of heavy leptons per unit dis- 
tance traveled in matter (in g/cm 2 ) is often approxi- 
mated as —(dE/dX) a + bE. The constant part, 
due to ionization losses, may be calculated using Bethe- 
Bloch formula, while the radiative part approximately 
proportional to E is due to a combination of stochastic 
e + e~ pair production, bremsstrahlung and photonuclear 
effects. In the EHE regime, the ionization loss term a is 
negligible compared to the stochastic losses. The radia- 
tive energy loss parameter b = b(E) varies slowly with 
energy for extremely high energy leptons, primarily due 
to an increase in the photonuclear energy loss rate b pn at 
very high energies. Although the radiative energy losses 
are in fact due to a series of discrete stochastic events, at 
high energies these interactions occur frequently enough 
that they can be considered quasi-continuous, increasing 
the overall brightness of the lepton track. 

While e + e~ pair production and bremsstrahlung are 
the dominant energy loss channels for fi above the TeV 
scale, for r bremsstrahlung is negligible and photonuclear 
effects dominate at EHE. Tau photonuclear losses are 
comparable to electromagnetic losses all the way down 
to the ionization-dominated region below E T ~ 20 TeV 

mm. 

Photonuclear energy losses by EHE leptons (and to 
a lesser extent bremsstrahlung losses) are not precisely 
known. Measurements of photonuclear cross sections 
from colliders must be extrapolated to very small x, 
and there are several models available in the literature 
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Different models of the nu- 
clear structure function are also available 0, m, l30| . 
The predicted loss rates for r leptons in ice are shown in 
Fig.^ Numerical values for energy loss rates given in this 
paper were evaluated using the MMC software package 
using ice as the default detection medium. 




Tau Energy [GeV] 

FIG. 1: Photonuclear energy loss per unit energy (6 pn ) for r 
leptons in ice, according to various models. BB/BS refers to 
|24|. plus the hard component of l2a|. Kok and ZEUS include 
the photon- nucleon cross sections of |26| and [2^|. respectively, 
instead of that from ALLM91 and ALLM97 refer to |H| 
and H3, with BM following [H Hj| but using the nuclear 
structure function from [3p|l . 

Recent calculations of muon bremsstrahlung [33l l34l 
l35l m| are in close agreement with each other (alth oug h 
they are about 5% higher than older calculations |37|h 
but disagree by as much as 20% for EHE taus. Because 
bremsstrahlung is suppressed by 1/mf, however, the dif- 
ferences between models have very little effect on the 
overall brightness of r tracks. 

B. Photonuclear Interactions 

In addition to the uncertainty in photonuclear energy 
loss rates at high energies noted above, there are several 
differences between photonuclear and electromagnetic in- 
teractions that must be considered. 

Both bremsstrahlung and pair production produce 
electromagnetic showers in the Cherenkov medium, with 
7's converting to e + e~ pairs and the e + and e~ in turn ra- 
diating more 7's. Photonuclear interactions, on the other 
hand, disrupt the nucleon involved and produce showers 
of hadrons which go on to interact with other nucleons 
in the medium. 

The light produced in a shower comes from the 
Cherenkov radiation of the many secondary particles in- 
volved in the shower. The total light yield is propor- 
tional to the integrated track length of the relativistic 
particles. The yield per unit shower energy is lower in 
hadronic showers because heavy particles have a higher 



threshold for Cherenkov radiation, energy is lost to the 
binding energies of the hadrons involved, and invisible 
slow neutrons are produced |3S| . 

The ratio of the light yield per unit energy in hadronic 
showers to that in electromagnetic showers depends on 
energy, as well. Neutral 7r mesons produced in the shower 
will decay to 77 and produce electromagnetic subshow- 
ers, but very few hadronic particles are formed in elec- 
tromagnetic showers. The production of 7r° is thus a 
"one-way street" |38l ] carrying energy out of the hadronic 
sector and into the electromagnetic. Roughly 30% of the 
energy in hadrons will go to 7r° in each generation of the 
shower, so high energy hadronic showers with more gen- 
erations will be more like electromagnetic showers. For 
small (~ 10 GeV) showers the ratio of light yields is about 
65%, rising to about 85% for a 10 TeV shower [H Hj 
and asymptotically approaching 100%. In the calcula- 
tions presented in Sect. ITTT1 we have assumed an average 
value of 75% for the energy deposited in hadronic pho- 
tonuclear interactions along the r (and fj,) track. 

Another important effect to consider is the energy 
spectrum of individual stochastic interactions. Relatively 
large fractions v of the total lepton energy can be de- 
posited in a single photonuclear interaction. In the anal- 
ysis of Ref. |25j in terms of a hard (perturbative) and a 
soft (non-perturbative) component, the energy spectrum 
v da/dv of the individual photonuclear interactions com- 
prising the perturbative component peaks at v ~ 10 -2 
for r leptons. Since the radiative energy loss parameter 
b < 10~ 6 g" 1 cm 2 , such interactions will be separated by 
several kilometers. 

Most often, such perturbative interactions will not 
occur along the contained track segment even in a 
kilometer-scale neutrino telescope. In this case, the per- 
turbative component of the photonuclear term (which 
in the analysis of Ref. [2j| dominates above ~ 10 PeV) 
should be ignored in estimating the brightness of the tau 
track, and only the non-perturbative component should 
be taken into account. If a perturbative interaction does 
occur within the detector, it will appear as a distinct 
cascade along the track, and it should still be possible to 
measure the "baseline" track brightness away from the 
shower (although the bright shower could confuse event 
analysis algorithms and thus lower the overall efficiency 
for detecting such events). 

As the tau energy increases over the EHE regime, how- 
ever, v da/dv shifts to lower v and the total photonuclear 
energy loss parameter b pn increases. Perturbative pho- 
tonuclear interactions thus become somewhat softer and 
much more common, with separations of a few hundred 
meters. In this regime, the quasi-continuous approxima- 
tion may be appropriate for the perturbative component 
as well as the non-perturbative part, although the details 
will depend on the particular detector and event recon- 
struction algorithm under consideration. In Sect. ITTT1 we 
present results for the two limiting cases, both of which 
should produce a detecatable signal: first assuming the 
perturbative photonuclear interactions of the r are very 
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rare, so that only the non-perturbative component of the 
photonuclear interactions contributes to the brightness 
of the r track; and then assuming the perturbative inter- 
actions are also quasi-continuous, increasing the average 
brightness of the r track and thus decreasing the ratio of 
the brightness of the tau to that of the muon. 



C. Radiative Decays 

Tau decays to /i are sometimes accompanied by initial- 
state or final-state radiation. This radiation follows a 1/fc 
spectrum at low energies (below the r mass scale), and 
is measured to occur with a branching fraction of 0.36% 
with a threshold of E y > 10 MeV in the r rest frame 
0. These photons will be boosted by a factor of up to 
7 (1 + 0), depending on the emission angle of the photon 
in the r rest frame, so that for E T = 1 PeV a photon 
with 10 MeV in the rest frame will have at most ~ 20 
TeV in the detector frame. Neutrino telescopes generally 
have an energy threshold of several TeV to a few tens of 
TeV for reconstruction of cascades (e.g. v e CC events), 
so 20 TeV is a reasonable benchmark for the energy at 
which a shower from initial-state or final-state radiation 
might be noticeable along a lepton track. For r tracks at 
the PeV scale, we therefore expect a sizeable shower at 
the decay vertex in 2% or fewer of observed tau decays. 
Because of the higher boost factor, noticeable showers 
will be about twice as common for EeV r decays. It is not 
clear that these showers would be distinguishable from 
the background of stochastic energy losses, but it might 
be helpful to treat sizeable showers along an observed 
track as candidate positions for tau decay vertices. 



III. DETECTABILITY 

An EHE r decaying to jj, will appear as a track which 
suddenly increases in brightness. The magnitude of the 
increase, taking into account the most probable fraction 
of the tau energy carried by the muon, the average en- 
ergy loss rates of the two leptons, and the relative light 
yields of hadronic and electromagnetic showers, is shown 
in Fig. |2 Losses to ionization, pair production and 
bremsstrahlung are also included; the bremsstrahlung 
model of [24| is used, but the results are nearly inde- 
pendent of the choice of bremsstrahlung model. 

Although the IceCube collaboration has not published 
any estimate of track energy resolution for the IceCube 
detector, a resolution of a (log 10 E^) ~ 0.3, correspond- 
ing to a factor of 2 in E^, was claimed for AMANDA-II 
[4l|. (In the radiative-dominated regime, the brightness 
scales approximately linearly with muon energy, so the 
resolution in brightness should be the same as the energy 
resolution.) One would expect IceCube to do at least this 
well, given the larger detector volume and better optical 
module electronics. 
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FIG. 2: The magnitude of the increase in brightness as a r 
decays to fj,, as a function of the r energy. The fi is assumed 
to take 1 — (y) — 0.4 of the r energy. The various lines 
correspond to different models for photonuclear energy loss, 
as in Fig.Q In the region of interest from 1 PeV to 1 EeV, the 
brightness steps up by a factor of between 3 and 7, depending 
on the model and on the energy of the r. 



The ANTARES collaboration expects a track ene rgy 
resolution of cr(log 10 E^) ~ 0.3-0.4 for PeV muons [42] 
(resolution at higher energies was not given). A review 
of the literature did not produce published track energy 
resolutions for the Baikal, NESTOR, or NEMO detec- 
tors. These smaller detectors (like AMANDA) do not 
have the effective volume to yield significant event rates 
at and above the PeV scale in any case, but the energy 
resolution for km 3 detectors in water should be compa- 
rable to, or better than, the ANTARES resolution. 

It appears that, for both ice and water km 3 neutrino 
telescopes, the expected track energy resolution should 
be sufficient to distinguish the brightness of the initial 
t track from that of the final fi track. We note that 
the energy resolutions given above refer to measuring 
a single energy for a through-going track, rather than 
trying to make separate energy measurements for dif- 
ferent segments of a track. However, in a km 3 detec- 
tor, the observed tracks will be several times longer than 
those visible in smaller instruments such as AMANDA or 
ANTARES, so enough information should be recorded for 
measurements of comparable accuracy. 

At the lower end of the EHE scale considered here (a 
few PeV) , the r track will only be visible if the neutrino 
interaction vertex is contained within the detector; even 
then, the r track will be quite short. Reconstruction of 
such events will be further complicated by the hadronic 
shower produced at the vN vertex, which will have typ- 
ically an energy of 0.25i?„. The ability to measure the 
energy of a short track segment within such a shower 
must be determined with a detailed and detector-specific 
Monte Carlo study, but it seems likely that a r track with 
a length of at least 200 m will be required. This would 
correspond to a neutrino energy threshold of around 5 
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ing only the final bang, known as a lollipop event, has 
already been pointed out. In this paper, we have shown 
that it should also be possible to tag taus which decay to 
muons. 

We note that this signature of r decaying to \i was con- 
sidered in Ref. |2^| , which concluded that the increase in 
brightness would not be detectable in a neutrino tele- 
scope. Two factors lead us to the opposite conclusion: 
the typical fraction of energy transferred to the fj, is some- 
what higher than they assumed, and the fact that a r 
loses energy primarily through photonuclear interactions 
reduces the brightness of the t track relative to that of a 
fi losing energy at the same rate. 



FIG. 3: The increase in brightness as r decays to \i, under the 
alternative assumptions that: (QC) the hard component of 
the photonuclear interaction is quasi-continuous, or (Non-C) 
that it is rare enough that it does not affect the measurement 
of the track brightness. The ability to resolve such showers 
and differentiate them from the track will depend on the neu- 
trino detector under consideration. 



PeV for detection of t using this signature, comparable 
to but slightly higher than the threshold for the double 
bang signature. The advantage of the muonic signature is 
that it should be remain detectable to very high energies, 
at least up to the EeV scale. 

As discussed in Sect. Ill 51 the curves in Fig. [21 assume 
that the light emission from photonuclear interactions 
is quasi-continuous, which may not be a valid assump- 
tion. Figure shows, for the model of [2j|, the effect of 
assuming to the contrary that hard photonuclear inter- 
actions are sufficiently rare that they do not contribute 
to the measured brightness of the underlying track. This 
assumption would hold if no such interactions occurred 
within the detector volume. Alternatively, an analysis 
searching for muonic tau decays could attempt to iden- 
tify bright showers along the observed track, and measure 
the "baseline" track brightness away from such showers; 
because very hard stochastic interactions are more com- 
mon for t than for fi, such an approach should heighten 
the contrast between fj, and r tracks. To our knowledge, 
no study of the ability of any neutrino telescope to re- 
solve showers along a track in this manner has yet been 
published. 



IV. DISCUSSION 



We also note the observation made in Ref. |2^| that the 
double bang and lollipop signatures are not completely 
free of experimental backgrounds, because muons which 
decay in flight can mimic these signatures. The authors 
of Ref. [23 estimate a rate of up to 50 km -3 yr _1 of 
these events. By contrast, there is no apparent physi- 
cal background to a track which suddenly increases in 
brightness as in the r — > \ivv decay signature. Exper- 
imental backgrounds will of course arise due to the in- 
trinsic variations of lepton energy deposition and photon 
detection which contribute to the detector energy reso- 
lutions quoted above. Because these variations depend 
on the particular detector and analysis techniques used, 
Monte Carlo studies will be needed to quantify the back- 
ground levels to be expected in actual experiments. 

At lower (TeV-PeV) energies than considered here 
thus far, r — > /1 might also be detectable if the neu- 
trino interaction vertex occurs within the detector, be- 
cause the two neutrinos produced in the r decay would 
carry off approximately half of the lepton energy. Al- 
though the r track would not be observed directly, the 
presence of the tau could be inferred on a statistical basis 
from the lower-than-expected energy of the muon track, 
when compared to the energy of the hadronic shower at 
the vN vertex. However, this signature would suffer from 
the same backgrounds as the 'inverted' lollipop (in which 
the tau production vertex and tau track are observed, 
rather than the tau track and decay vertex). These sig- 
natures can be faked by a v^N CC interaction with high 
y, where less than the mean energy is transferred to the 
outgoing /j,, or by CC v e or NC v x interactions where 
secondary lower-energy /i are produced via ir^ decay in 
the hadronic shower. 



The detection of tau neutrinos offers an excellent 
method for detecting astrophysical neutrinos, since there 
is essentially no atmospheric tau neutrino background. 
However, the classic double bang signature of tau neutri- 
nos is only expected to be detectable over a single decade 
of energy, from a few PeV to perhaps 20 PeV. The possi- 
bility of tagging tau events at higher energies by detect- 



We believe the signature of muonic tau decay will be 
useful in identifying astrophysical tau neutrino events in 
the coming generation of kilometer-scale Cherenkov neu- 
trino telescopes such as IceCube. This signature may be 
particularly important in the energy region above a few 
tens of PeV, where the classic double-bang signature is 
no longer observable. 
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